Structural study of pyrimidine nucleoside analogues. I: Molecular mechanics and semiempirical calculations of 2′-deoxy-2′-fluoroarabinofuranosyluracils

Molecular mechanics (MM) calculations have been performed on the title compounds. For the MM minimum energy conformation obtained by conformational analysis, molecular orbital (MO) calculations (MNDO and AM1) have also been performed. The geometries obtained have been compared with the experimental ones extracted from the Cambridge Structural Database (CSD). A qualitative structure-activity relationship has been pointed out based on the electrostatic potentials calculated at different positions on the electronic surface.     

Role of the tautomerism of 2-azaadenine and 2-azahypoxanthine in substrate recognition by xanthine oxidase

The tautomerism of 2-azaadenine and 2-hypoxanthine has been examined in the gas phase and in aqueous solution. The tautomerism in the gas phase has been studied by means of semiempirical and ab initio quantum-mechanical computations, as well as density-functional calculations. The influence of the aqueous solvent on the relative stability between tautomers has been estimated from self-consistent reaction field calculations performed with different high-level continuum models. The results provide a detailed picture of the tautomeric preference for these purine bases. The importance of tautomerism in the substrate recognition by xanthine oxidase is discussed. Finally, the rate of oxidation of 2-azaadenine and 2- hypoxanthine by xanthine oxidase is discussed in terms of the recognition model at the enzyme active site.     

Intraframework potential energy function of zeolites

Potential energy functions suitable for the (Cs_1–xNa_xT₂O₄)-A type zeolite family were obtained from the known crystal structure of Cs₇Na₅-A zeolite. Using these potential functions, several other crystal structures in these family were obtained by molecular mechanical calculations. The activation energy - to -cage transmission of the H₂ molecule and the void volume of Cs₃Na₉-A zeolite were calculated.     

取代苯体系的二阶非线性光学性质:动力学李代数方法

提出了一种动力学李代数方法来研究取代苯体系的非线性光学性质. 对于给定的PPP模型(Pariser-Parr-Pople)哈密顿量, 生成了一个动力学李代数. 依据这些代数元构造出演化算子作为群参数的函数, 通过求解一组非线性微分方程能够得到这些群参数. 再按照统计力学中的密度算子公式给出取代苯分子体系偶极矩的统计平均值. 于是导出二阶极化率的表达式. 与其他量子力学计算结果比较, 表明这种动力学李代数方法在预言有机共轭分子的非线性光学性质上同样有用.     

表皮生长因子受体和抑制剂之间分子对接的研究

研究了表皮生长因子受体(EGFR)和4-苯胺喹唑啉类抑制剂之间的相互作用模式，表皮生长因子受体的三维结构通过同源蛋白模建的方法得到，而抑制剂和靶酶结合复合物结构则通过分子力学和分子动力学结合的方法计算得到。从模拟结果得到的抑制剂和靶酶之间的相互作用模式表明范德华相互作用、疏水相互作用以及氢键相互作用对抑制剂的活性都有重要的影响，抑制剂的苯胺部分位于活性口袋的底部，能够与受体残基的非极性侧链产生很强的范德华和疏水相互作用，抑制剂双环上的取代基团也能和活性口袋外部的部分残基形成一定的范德华和疏水性相互作用，而抑制剂喹唑啉环上的氮原子能和周围的残基形成较强的氢键相互作用，对抑制剂的活性有较大的影响，计算得到抑制剂和靶酶之间的非键相互作用能以及抑制剂和靶酶之间的相互作用信息能够很好地解释抑制剂活性和结构的关系，为全新抑制剂的设计提供了重要的结构信息。     

Structure and torsional flexibility of the linkage between guanine and fluorene residues in the deoxyguanosine–aminofluorene and deoxyguanosine–acetylaminofluorene carcinogenic adducts

Potential energy surfaces for rotations around two central CN bonds in N-(deoxyguanosin-8-yl)-2-acetylaminofluorene (AAF–dG) and its deacetylated derivative (AF–dG) were studied using Amber 95 molecular mechanics. Both of these adducts are known to be strong mutagens and carcinogens. New Amber 95 force field parameters were derived for the linkage connecting guanine and fluorene moieties in AAF–dG and AF–dG. For this purpose, we determined ab initio MP2/cc-pVDZ//B3-LYP/6-31G* and polarized continuum model Hartree–Fock/6-31G* potential energy surfaces of smaller model systems that included the N-methylimidazole–acetylaniline and N-methylimidazole–aniline adducts. The molecular mechanics parameters were adjusted to minimize differences between the gas-phase ab initio and molecular mechanics surfaces of these model systems. The resulting parameters were transferred to AF–dG and AAF–dG. The barrier for the rotation of the fluorene residue in AF–dG was found to be less than 2 kcal/mol. Such a small barrier renders the fluorene moiety freely rotatable at room temperature. In contrast, the fluorene rotation in AAF–dG is hindered by a significantly larger barrier of 10 kcal/mol. This barrier corresponds to conformations in which the fluorene and acetyl groups lie in the same plane, and is largely due to steric repulsion. Similarly, the coplanar arrangement of guanine and the bridging amino or acetyl groups is disfavored by 5–10 kcal/mol, with AAF–dG again being the more rigid of the two molecules. Energy minima for a rotation around a bond between guanine and the bridging nitrogen are found at ±80° in AAF–dG, and at 120° and –90° for AF–dG. Overall, the fluorene–dG linkages in AF–dG and AAF–dG adducts have significantly different equilibrium structures and torsional flexibilities. These differences may be contributing factors for the observed disparity in mutagenic effects of these adducts.Electronic Supplementary Material: Supplementary material is available in the online version of this article at Acknowledgements. This work was supported by the NSF REU grant no. CHE-0243825 to Loyola University Chicago. We thank to Tom Ellenberger and Shuchismita Dutta for providing us with their results prior to publication.     

Chromium(III) Hydroxide Solubility in the Aqueous K+-H+-OH?-CO2-HCO 3? -CO 32? -H2O System: A?Thermodynamic Model

Chromium(III)-carbonate reactions are expected to be important in managing high-level radioactive wastes. Extensive studies on the solubility of amorphous Cr(III) hydroxide solid in a wide range of pH (3–13) at two different fixed partial pressures of CO₂(g) (0.003 or 0.03 atm.), and as functions of K₂CO₃ concentrations (0.01 to 5.8 mol⋅kg⁻¹) in the presence of 0.01 mol⋅dm⁻³ KOH and KHCO₃ concentrations (0.001 to 0.826 mol⋅kg⁻¹) at room temperature (22±2 °C) were carried out to obtain reliable thermodynamic data for important Cr(III)-carbonate reactions. A combination of techniques (XRD, XANES, EXAFS, UV-Vis-NIR spectroscopy, thermodynamic analyses of solubility data, and quantum mechanical calculations) was used to characterize the solid and aqueous species. The Pitzer ion-interaction approach was used to interpret the solubility data. Only two aqueous species [Cr(OH)(CO₃)₂²⁻ and Cr(OH)₄CO₃³⁻] are required to explain Cr(III)-carbonate reactions in a wide range of pH, CO₂(g) partial pressures, and bicarbonate and carbonate concentrations. Calculations based on density functional theory support the existence of these species. The log ₁₀ K° values of reactions involving these species [{Cr(OH)₃(am) + 2CO₂(g)⇌Cr(OH)(CO₃)₂²⁻+2H⁺} and {Cr(OH)₃(am) + OH⁻+CO₃²⁻ ⇌Cr(OH)₄CO₃³⁻}] were found to be −(19.07±0.41) and −(4.19±0.19), respectively. No other data on any Cr(III)-carbonato complexes are available for comparisons.     

Multiple glass transitions and local composition effects on polymer solvent mixtures

Recent differential scanning calorimetry (DSC) results on polystyrene–solvent mixtures show two distinct glass transitions whose positions and widths vary with composition. Parallel work on the dynamic response in polymer blends has focused on how segmental mobilities are controlled by local composition variations within a “cooperative volume” containing the segment. Such variations arise from both chain connectivity and composition fluctuations. We account for both using a lattice model for polymer–solvent mixtures that yields the composition distribution around polymer and solvent segments. Insights from our lattice model lead us to simplified calculations of the mean and variance of local composition, both in good agreement with lattice results. Applying our model to compute DSC traces leads to an estimate of the cooperative volume, since a larger cooperative volume both reduces the biasing effect of connectivity, and narrows the composition distribution. Comparing our results to data, we are able to account for the composition-dependent broadening with a cooperative length scale of about 2.5 nm. © 2006 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys 44: 3528–3545, 2006     

Superstability estimates for anharmonic systems

We consider an anharmonic crystal described by variablesS _x ,x ^d ,S _x , with one-body interaction ¦S _x ¦ and nearest neighbor (n.n.) two body interaction ¦S _x –S _y ¦. We prove that, for ^d bounded, , where is the correlation function for the free boundary condition Gibbs state in ,>0 and are suitable constants independent of and . This generalizes previous results obtained in the case.Research partially supported by Consiglio Nazionale delle Ricerche.     

An evaluation of molecular models of the cytochrome P450 Streptomyces griseolus enzymes P450SU1 and P450SU2

Summary P450SU1 and P450SU2 are herbicide-inducible bacterial cytochrome P450 enzymes from Streptomyces griseolus. They have two of the highest sequence identities to camphor hydroxylase (P450cam from Pseudomonas putida), the cytochrome P450 with the first known crystal structure. We have built several models of these two proteins to investigate the variability in the structures that can occur from using different modeling protocols. We looked at variability due to alignment methods, backbone loop conformations and refinement methods. We have constructed two models for each protein using two alignment algorithms, and then an additional model using an identical alignment but different loop conformations for both buried and surface loops. The alignments used to build the models were created using the Needleman-Wunsch method, adapted for multiple sequences, and a manual method that utilized both a dotmatrix search matrix and the Needleman-Wunsch method. After constructing the initial models, several energy minimization methods were used to explore the variability in the final models caused by the choice of minimization techniques. Features of cytochrome P450cam and the cytochrome P450 superfamily, such as the ferredoxin binding site, the heme binding site and the substrate binding site were used to evaluate the validity of the models. Although the final structures were very similar between the models with different alignments, active-site residues were found to be dependent on the conformations of buried loops and early stages of energy minimization. We show which regions of the active site are the most dependent on the particular methods used, and which parts of the structures seem to be independent of the methods.